Classical decays in decoherent quantum maps 
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We study the asymptotic long-time behavior of open quantum maps and relate the decays to the 
eigenvalues of a coarse-grained superoperator. In specific ranges of coarse graining, and for chaotic 
maps, these decay rates are given by the Ruelle-Pollicott resonances of the classical map. 
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The study of the emergence of classical features in the 
quantum behavior of hamiltonian systems has centered 
on two areas: on the one hand, the statistical fluctuation 
properties of the quantum spectrum were the first quan- 
tities to reveal universality classes related to classical - 
chaotic or integrable - behavior On the other hand, 
the study of time dependent quantities has also provided 
quantum signatures of the classical world. The growth 
of Von Neuman entropy in decohering quantum systems 
0, Q showed striking differences separating chaotic and 
regular systems and pointed to the relevance of Lyapunov 
exponents - more precisely the Kolmogorov-Sinai entropy 
- for this quantity's time behavior. Recently there has 
been a strong interest in the study of the Loschmidt echo 
Q - specially in the context of fidelity decay in quantum 
maps and algorithms - which displayed essentially the 
same characteristics. In this letter we report the emer- 
gence and the relevance in quantum decohering maps of 
yet other quantities of a classical nature - the Ruelle- 
Pollicott resonances 0, @ - and show that they rule the 
asymptotic decay of time dependent quantities. 

In [2j a simple model for decoherence in quantum maps 
was introduced based on the Kraus representation of su- 
peroperators ||, with a straightforward interpretation 
as quantum coarse-graining in phase space. The model 
was employed to study the entropy evolution and to con- 
firm the prediction [2j that for certain regimes the rate 
of growth is independent of the coarse graining and is 
given by the Lyapunov exponent of the classical motion. 
This is a short time regime, of the order of the Ehrenfest 
time, after which relaxation towards the uniform density 
occurs. In this letter we study this asymptotic relaxation 
regime for quantum maps and show that it is ruled by 
the eigenvalues of a similar coarse-graining superopera- 
tor. These eigenvalues can be effectively calculated and 
they determine - in certain well defined ranges - classical 
decay rates independent of the coarse graining. These 
decay rates are given, in the semiclassical limit of chaotic 
systems, by the Ruelle-Pollicott resonances. We show in 
particular that, beyond the Ehrenfest time, these reso- 
nances determine the asymptotic behavior of the linear 
entropy and of the Loschmidt echo. 

The unitary evolution of the density matrix for a quan- 
tum system is given by p n +i = U p n U^ where U is a uni- 
tary evolution operator that quantizes a classical map. 



For simplicity we consider area preserving maps on a 
phase space with periodic boundary conditions, i.e. a 
torus (normalized to unit area). In this case the Hilbcrt 
space has a finite dimension N related to h by 2irhN = 1 . 
Position and momentum eigenstates are respectively: 
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i,j€[0,N-l]. 
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They form an N x N grid that constitutes the quantum 
phase space for these maps and their overlap is the dis- 
crete Fourier transform (DFT) of dimension TV. Transla- 
tions on this grid can be implemented by unitary opera- 
tors Tip, q) , easily constructed from cyclic shifts and the 
DFT 9]. They carry the integer labels p, q and they form 
a complete unitary operator basis in Hilbert space. The 
linear action of U on the vector space of density matrices 
defines a superoperator L = U <£> W , such that 
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which is still unitary, of dimension N 2 x iV 2 . The eigen- 
values of L are e^ Ei ~ Ej \ where Ei are the Floquet eigen- 
values of the map. 

The quantum coarse graining is implemented by a 
Kraus superoperator 
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De=Y. c e ip,q)f(p,q)®Pip,q), 



(3) 



p,q=0 



where c e (p, q) is a smooth positive periodic function on 
the torus to be defined below. In order for D e to be trace 
preserving the condition c e (p, q) = 1 must be satis- 
fied. Inasmuch as c e ip,q) is a narrow gaussian-like func- 
tion, D e implements an incoherent sum of slightly dis- 
placed density matrices over a phase space region of order 
e, this being a measure of the amount of coarse graining. 
To avoid an overall drift we further assume that c £ (p, q) is 
an even function of its arguments. The spectral proper- 
ties of D e are easily calculated. D e is hermitian and using 
the fact that f(p,q)f(\,n)fi(p,q) = T(A, /^e^C^-m) 
its spectral properties are given by 



D e f(\,n) =c-(A,m)T(A, M ), 



(4) 



where c~(A,^) = T, p , q Ce{p,q)e i ^' iXq ^ p) is the DFT of 
c t ip,q). It is irrelevant whether we specify the diffusion 
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by c e (p,q) or c e (A,/i). We found more convenient to de- 
fine the latter as 



c e (A,/i) 



e - 1 ( ^ ) 2 (sin 2 [7rA/AT]+sin 2 [tt^/AT] ) 



(5) 



whose Fourier transform is a smooth periodic gaussian- 
like coarse graining kernel, of width e/27r. The coarse- 
grained dynamics of the quantum map will be given by: 



Pn+l = O e Lp n = L e p n . 



(6) 



For any finite value of e , L e is a convex sum of uni- 
tary matrices and is therefore a contracting map whose 
spectrum is contained in the unit circle, and has a non 
degenerate eigenvalue equal to 1, corresponding to the 
uniform eigenvector T(0,0), and (iV 2 — 1) eigenvalues of 
modulus smaller than 1. Therefore, the time evolution of 
the density matrix can be decomposed as a sum of expo- 
nentially decaying modes plus a constant. In particular, 
the asymptotic decay towards the uniform density will 
be given by the eigenvalue closer to the unit circle. Al- 
though this is strictly true for any value of N and e > 0, 
the behavior of the eigenvalues in the limits N — > oo 
and e — ► is strongly dependent on the classical features 
-chaotic or regular - of the map in question. 

Before studying these limits we review the equivalent 
procedure for classical maps: densities in phase space 
evolve with the Frobenius-Perron (F-P) operator £ aris- 
ing from the Liouville equation. For a classical map, 
the general form of C is £(x, y) = S[y — /(x)] where 
x = (q,p), y = (q',p')- As it stands, the operator C is 
unitary on L 2 , the space of square integrable functions 
on phase space. When it is convoluted with a narrow 
coarse graining operator T> ei e being the coarse graining 
parameter, we obtain 



£ £ = V e o L. 



(7) 



T> e has the role of damping the high frequency compo- 
nents of phase space distributions in L 2 thus produc- 
ing an effective truncation of C and drastically chang- 
ing its spectrum. In fact now C e is compact and its 
spectrum consists of isolated resonances contained in the 
unit circlefTol [T]|. The behavior of these resonances as 
e — > is radically different for chaotic or regular maps. 
If the map is ergodic and mixing, C e has an isolated 
non-degenerate eigenvalue equal to 1, corresponding to 
the uniform density, and isolated eigenvalues with finite 
multiplicity, and modulus strictly smaller than 1. In the 
singular limit e — > 0, these eigenvalues are the Ruelle 
resonances |0, |n| and they determine the exponential 
decay of correlation functions in phase space [f| El • An 
equivalent procedure to uncover these resonances consists 
in introducing a basis of functions ordered by resolution 
in phase space, then constructing C in this basis and sub- 
sequently truncating to a finite dimension M, resulting 



in a matrix C , whose larges t eig envalues stay fixed as 
M -> oo (see El E! El El, E3 for a review and E3 
for a quantum application). 

It was proved recently, under quite broad conditions 
E3 (Theorem I), that for smooth maps on the torus 
and for any fixed e the spectrum of the coarse grained 
quantum propagator L e converges to that of C e in the 
semiclassical limit N — > oo . This in turn implies that at 
finite e the quantum decay rates are given by the corre- 
sponding classical eigenvalues. We study below the issue 
of letting e — > 0. 

We first compare these spectra for the perturbed 
Arnold cat map (and its quantization |18() 



P' = P + <1 - 27rfcsin[27T(?] , , 
q'= q + p / + 2nksm[2np'}. { a >' 



(8) 



The classical Arnold cat map is ergodic and mixing, has 
uniform hyperbolicity and its Lyapuonov coefficient is 
X = (1/2) (3 + y/E). The perturbation is added to avoid 
the non-generic behavior of the quantized version and \ is 
not changed substantially by it. On the contrary, the Ru- 
elle resonances are very sensitive to the perturbation. In 
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FIG. 1: Comparison of the leading eigenvalues of the quan- 
tum L e and classical C e coarse-grained propagators for the 
perturbed Arnold cat map, with different coarse graining and 
perturbation (a) k = 0.01, e = 0.4; (b) k = 0.02, e = 0.314. 
(• = classical, grid of 50 x 50 ; ffi = quantum, TV = 150). 



Fig^we compare the eigenvalues of C e and L e . The main 
difficulty in both calculations is the large dimensions of 
the matrices involved. For the classical calculation we 
took the coarse graining operator T> e to be a periodic 
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gaussian of width c/(2tt). The spectrum of £ £ was ob- 
tained by projecting eq. © on a grid of 50 x 50 sites and 
diagonalizing. The procedure is stable for the leading 
eigenvalues as long as e is much larger than the grid size 
[lOj. The spectrum of L e was obtained by a variational 
method similar to that proposed in 0] . A basis of right 
and left trial functions, which are smooth on the unstable 
and stable manifold respectively, are used to construct a 
subspace that contains the leading eigenspace. The ad- 
vantage of this method is that the resulting generalized 
eigenvalue problem has a small dimension as the basis 
contains information about the dynamics. Details about 
the computation of the spectrum of L, with this method 
will be reported in a future article [20j • Good agreement 
for the leading eigenvalues of both spectra was found in 
all the cases we computed, provided N was large enough, 
in agreement with . In order to obtain the Ruelle res- 
onances we now have to take the limit e — * 0. Clearly 
this limit cannot be taken at constant N because for finite 
matrices unitarity is eventually recovered and all eigen- 
values go back to the unit circle. 

In Fig. |2 we study this process as a function of N and 
e for the leading eigenvalue Ai. We observe that for each 
N there is a range of values of e for which the eigenvalue 
is independent of e and for which the limit e — > can be 
extrapolated safely, thus defining a property of the classi- 
cal map, independent of N and e. A similar situation oc- 
curs for the other eigenvalues, but the safe range becomes 
smaller as the distance to the unit circle increases. We 
conclude therefore that the leading portion of the spec- 
trum of L e for large N is of a purely classical nature and 
provides the rates of asymptotic decay of quantum time 
dependent quantities. For chaotic maps these classical 
rates are further identified with the spectrum of Ruelle 
resonances. |21| 
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FIG. 2: The leading eigenvalue |Ai| as a function of e, for dif- 
ferent values of N. There is a range, growing with N where 
the eigenvalue is independent of e so that it can be extrapo- 
lated to e -y 0. k = 0.01. 



To test this prediction consider first the evolution of 
the linear entropy S n = — lntr[p„] . It was shown in 
0, that for short times the entropy growth is linear 
with a slope independent of e and determined by the 
Lyapunov exponent of the map. This is shown in the 
inset of Fig. [21 where S n is plotted as a function of n for 
N = 450, e = 0.05 and various values of the perturbation. 
As expected for an ergodic map, for longer times the 
state relaxes to the uniform distribution and S converges 
to Soo — — ln(l/AT). To uncover the rate of relaxation 
we subtract the uniform density from the initial state 
and obtain, using the long time decomposition of p n into 
decaying modes, S n — —2 In | Ai \n. Again a linear growth 
is obtained but this time with a slope dependent on the 
leading Ruelle resonance. This feature is clearly seen in 
Fig. |3| where we plot S n for different values of k. 
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FIG. 3: Linear entropy S n for the perturbed Arnold cat for 
several values of k. The initial value is a coherent pure state 
with the uniform density subtracted. N = 450, e = 0.05. No- 
tice the two linear regimes: the short time Lyapunov (inset) 
and the asymptotic Ruelle. Solid lines are computed from the 
leading eigenvalue of L e . 

A very similar pattern emerges for the decay of the 
so called Loschmidt echo, or fidelity This is de- 

fined as M n — tr(p n p' n ) where p n and p' n are states 
evolved with slightly different maps from an initial state 
po- This quantity, originally considered by Peres as 
a measure of sensitivity to perturbations, has aquired 
renewed importance as a measure of environment inde- 
pendent decay of spin echoes and of the accuracy 
of implementation of quantum algorithms. It is usually 
considered only in the case of strictly unitary evolution 
where it equals M„(0) = tr[L"po L''Vo]- For (slightly) 
open systems it is natural to replace this quantity by 
M„(e) = tr[L"p p ] and then consider the limit 
e — > to extract features that are independent of the 
coarse-graining. The same considerations regarding the 
non-commutativity of the N and e limits apply here and 
in the "safe" regime of N and e the asymptotic decay 
of M n is again ruled by the Ruelle resonance. In the 



4 



classical similar result for the echo was obtained 

recently j^. Our discussion shows that the quantum 
decay must follow the classical one. In fact the echo 
is closely related to the linear entropy discussed above. 
Schwartz inequality tr(p n p' n ) < y/ tr(/?2 )tr(p'^) leads to 
lnM„ < -(l/2)(5 n + S' n ). Again the uniform density 
is subtracted from the initial state to reveal the expo- 
nential relaxation regime. In the asymptotic regime the 
inequality becomes an equality and we find that beyond 
the Ehrenfest time the decay proceeds as 

lnM„ = -(l/2)(ln|Ai| +ln|Ai|)n. (9) 

This is illustrated in Fig. where it is clearly seen that, 
after the initial Lyapunov regime the echo decays with 
the average of the two linear entropies. 
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FIG. 4: Decay of the Loschmidt echo beyond the Ehrenfest 
time. Two values of &,(+ = k = 0.005, * = k = 0.017 per- 
turbed by Afc = 0.002 are plotted. N = 450, e = 0.1. The full 
line shows the echo while the dotted lines give the perturbed 
and unperturbed entropies — S n and — S' n . 

Our spectral analysis shows that for any observable 
whose quantum evolution depends on the operator L, and 
this is not restricted to maps only, the best procedure to 
extract its classical behavior is to replace L by L e and let 
e — ► and N — > oo (meaning h — > ) in such a way as to 
remain in the "safe" region where unitarity is not recov- 
ered. The convergence of the spectra of Land L e , proven 
under very specific conditions in 0, 0], and demon- 
strated here in numerical experiments, is expected to be 
much more general, and in practice provides a very pow- 
erful means of extracting classical behavior from quan- 
tum maps 25] . Whether different models of decoherence, 
based on physical mechanisms in small quantum devices, 
lead to the same results remains to be explored but, at 
least from the numerical point of view, are also subject 
to this spectral analysis and can reveal other modes of 
decay not related to classical properties. 
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